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Abstract: 


The recent acquisition of data on the distribution of galactic 
rray emission has made possible a reconsideration of the long-standing 
controversy concerning the existence and extent of the galactic halo. 
Analysis of the implications of the SAS-2 y-ray d’ta, making use of 
recent CO line emission and other data for determining the large-scale 
distribution of galactic gas, implies that there is a nonuniform distri- 
bution of cosmic rays in the galaxy. This fact rules out large trapping 
halo models and particularly the recently proposed closed halo models 
in the same way that it rules out extragalactic origin models. 

We also consider detailed models of diffusion halos of various 
sizes perpendicular to the galactic plane. In such models, the scale 
perpendicular to the plane has a strong effect in determining the radial 
distribution of cosmic rays. Such radial distributions are calculated 
for cylindrical coordinate models. The implied y-ray longitude distri- 
butions are then calculated and compared with the SAS-2 data for 
goodness-of-fit. Assuming the sources to be supernova remnants or pulsars, 
cosmic ray nucleon halo models with scale heights greater than 3 kpc 
are found to provide a poor fit to the y-ray longitude data (probability 
of 6 % or less). Thin halo, or sourct dominated diffusion models are 
found to provide a good fit to the y-ray data, with an upper limit scale 
height of 'v. 3 kpc. 

Consideration of the y-ray latitude data gives a half thickness 
of 2 + 2 kpc for the cosmic ray electron halo. A best-fit value of 
the product the galactic center, based on an estimate of 

the Compton production in the central region gives = 0.5 kpc. 


I Introduction: 


The concept of a large quasispherical region of galactic 
cosmic rays surrounding our Galaxy, i.e. a "halo" or "corona," 
has become the subject of much controversy in the last decade. 

Most of the discussion has centered upon the question of the 
existence of an extensive halo of galactic radio emission with 
some discussion centering upon the question of the effective 
"lifetime" or trapping time of cosmic rays in the galaxy. 

Shklovskii (1952) showed that the observed high latitude 
distribution of nonthermal radio emission could be naturally 
explained by postulating an extensive halo of radio emissivity 
surrounding the galaxy. Pikel'ner (1953) suggested that there 
should exist an accompanying gaseous halo and magnetic field of 
'V- 3x10*^ G in order to retain the cosmic rays producing the 
radio emission. Observations by Baldwin (1955) appeared to 
confirm the existence of an extensive radio halo, this halo 
having a radial dimension of the order of 20 kpc (Baldwin 1967). 

A somewhat flattened halo model with an axial ratio of '^'1.5 was 
suggested by Mills (1959). Biermann and Davis (1958) pointed out 
that an extensive cosmic ray halo would naturally account for the 
observed isotropy of cosmic rays. An electron intensity for the 
halo similar to that in the galactic disk was deduced by Felten (1966). 
However, the existence of the radio halo was already being questioned 
by Burke (1967). Some of the early radio work may have been affected by 
sidelobe contamination and the additional effects of radio spurs 
and loops (nearby features), which became apparent with higher 
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resolution. An extragalactic background also complicates the analysis. 
More recent analyses have questioned the existence of an extensive 1 
halo (Wielibinski and Peterson 1968, Yates 1968, Razin 1971, Ilovaisky 
and Lequeux 1972, Price 1974), or favored a weak non-confinino halo 
(Webster 1975) or a very thin halo or thick disk of a few kpc 
extent (Illovaisky and Lequeux 1972, Baldwin 1976). The opposite 
point of view in support of an extensive radio halo has been 
taken by Daniel and Stephens (1975 and references therein) and Bulanov 
Syravatskii and Dogiel (1976) and Ginzburg and Ptuskin (1976 
and references therein, see also Ginzburg and Syrovatskii 1964). 

There have also been searches for radio halos associated with 
edge-on spiral galaxies. Baldwin and Pooley (1973, see also Van 
derKruit and Allen 1976) have found that a flat halo (or flat disk) 
of half thickness 'V' 2 kpc around NGC 891. Ekers and Sancisi (1977) 
have found a more extensive but still flattened halo around 
NGC 4631, however, these authors point out that NGC 4631 may not 
be a good analogue of our Galaxy since (1) it has an unusually bright 
central emission region and (2) it is surrounded by large 
concentrations of HI and (3) it has a nearby disturbed companion 
galaxy. Searches for halos in other edge-on spirals have so far 
given no positive results (e.g. Mills 1957) and the question of the 
existence of a radio halo around M31 has been clouded in interpretational 
difficulties (van der Kruit and Allen 1976). Further compli- 
cating the situation is the fact that the lack of an extensive 
radio halo is not a conclusive argument against the existence 
of a large halo in which cosmic rays are confined (Ginzburg and 
Ptuskin 1976). 
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Measurements of the ratio of Be/B and the fraction of ^*^Be in 

the cosmic radiation have been made in order to determine the age 

of the cosmic rays from the survival uf ^^Be produced by cosmic 

ray spallation (Webber et al 1973, O'Dell et al 1973, Hagen, Fisher ar:d 

Ormes 1977, Garcia-Munoz 1975). Studies of cosmic ray secondaries (spalled 

nuclei and positrons) indicate that cosmic rays travel on the 

average through about 5 g/cm of matter (Shapiro, Silberberg and Tsao 

1975) A cosmic ray age of x 10^ yr would imply that the cosmic 

rays have traveled through a medium of average atomic density > 0.4 

cm and therefore favor a disc containment model or one in which cosmic 

rays diffuse rapidly in the halo so that once they leave the disk 

have little chance of returning in a diffusion model. On the 

other hand a determination of r >> 10^ yr would favor containment in 

-2 -3 

an extensive halo (n 10 cm ). The measurements quoted above 

appear as yet to provide an inconclusive or borderline test as the 

range of overlap at present appears to be (1-2) x 10^ yr. The 

data of the Chicago group (Mason 1977) at energies of 100 MeV/nucleon 

7 10 

give an age of 1.7 x 10 yr from Be measurements while the age 
in the range of several GeV/nucleon (more typical of "cosmic rays" 
and also the range which is important for TT-decay y-ray production 
(Stecker 1973)) is given by O'Dell et al . (1975) to have a 1 o 
upper 1 imit of 10^ yr. 

In view of the situation with the radio data, and the additional 
present uncertainty regarding the age of galactic cosmic rays (O'Dell 
et al. 1973, Hagen et al 1975, Garcia-Munoz et al . 1975), it would seem 
desirable to have an independent test of the existence and extent of 
a cosmic ray halo. Such a test is furnished by the recent observations 
of the distribution of galactic 100 MeV y-rays and their interpretation 
(Stecker 1976a, 1977). 3 


II. Definition of a Halo 

Subsequent to their initial suggestion of a galactic halo, Pikel'ner 

and Shklovskii (1957, 1958, 1959) developed a model for the halo which 

related the spherical system of older (population II) stars, a gaseous 

halo, a cosmic ray halo and the radio halo. The halo of population ^ 

II stars certainly exists. However, for the purposes of detailed 

consideration of the observational consequences and in view of the 

difficulties in interpretation, it may perhaps be better to separate 

out these various types of halo. Energy losses by synchrotron 

and Compton radiation deplete the cosmic ray electron component over 
7 8 

10 -10 yr and thus shrink the observable radio halo while still 
allowing for an extensive cosmic ray nucleon halo (see, e.g. 

Bulanov, Syrovatskii and Dogiel 1976, Ginzburg and Ptuskin 1976 
and references therein). Thus, evidence of an extensive radio 
halo may be taken as strong support for an extensive nucleon halo 
bit the absence of such evidence may not with certainty rule out 
the existence of a nucleon halo. Such evidence makes cosmic ray 
trapping unlikely (Webster 1975) but would not rule out diffusive 
type models. We are faced then with a situation of uncertainty 
with regard to the putative nucleonic cosmic ray halo. It is a 
discussion of this nucleonic halo (hereafter to be called the 
"cosmic ray halo" or "halo") to which to restrict ourselves by 
addressing the y-rhy data. Since cosmic-ray nucleons do not 
suffer significant energy losses in a halo of fairly extensive size, 
evidence against an extensive cosmir-ray halo also argues against 
an extensive radio halo. 
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III. The Uniform or Closed Halo Models 

Two types of halo models are usually considered, viz. 

(a) the uniform or closed type in which cosmic rays are trapped 
for long periods in the halo and can reenter the disk many times 
and (b) the diffusion halo. The closed halo model has been 
recently revived by Rasmussen and Peters (1976). The diffusion 
halo models have been extensively examined by various theorists 
at the Lebedev Institute (Bulanov, Syrovatskii and Dogiel, 1976, 

Ginzburg and Ptuskin 1976 and references therein). In addition, 
the uniform and diffusion halo models appear to have been confused 
in the literature in the past as has been pointed out by Ginzburg 
and Ptuskin (1976). The reader should keep the difference in mind 
in the subsequent discussion. 

It has been previously noted that the SAS-2 observations of 
100 MeV Y-radiation in the galaxy imply a nonuniform cosmic- 
ray distribution in the galaxy which argues against extragalactic 
cosmic ray origin (Stecker 1975, Dodds et al 1975, Stecker et al 1975). 

The implied existence of a nonuniform galactic cosmic-ray distribution 

puts new restrictions on halo models for propagation and containment 

of cosmic rays. If cosmic rays are produced primarily as a 

result of supernova explosions or in pulsars, the striking correlation 

of the galactic cosmic ray distribution and the supernova remnant 

distribution (Stecker 1975a) and pulsar distribution as given by 

Seiradakis (1976) implies source dominated diffusion (or convection, 

see Jokipii 1976) of cosmic rays out of the galaxy. This argues 

against trapping of cosmic rays in a large halo with multiple 

reflections of cosmic rays at the boundaries. Such trappings would produce a^ 




more uniform distribution of cosmic rays in the galaxy than the y-ray 
• observations imply (Stecker 1976a, 1977). The large halo model with 

O 

a long 10 yr) trapping time has been advocated in the past (see, 
e.g. Ginzburg and Syrovatskii 1964). An attractive feature of the large 
halo model has been that it eases the requirements for isotropization of 
cosmic rays in the solar galactic neighborhood without invoking extra- 
galactic cosmic ray origin (Biermann and Davis 1958). Recent analysis 
of the high galactic latitude radio data indicates however, that the magnetic 
field in the halo region may not be strong enough to trap cosmic rays 
(Webster 1975). 

The existence of any type of closed galaxy model with long term 

Q 

(> 10 yr) confinement of cosmic rays is brought into question by the 

Y-ray results. This includes the closed galaxy model (Rasmussen and 

Peters 1976) where cosmic rays have undergone spallation through 
2 

X = pet 'V- 100 g/ern . The isotopic measurements also argue against the 
uniform confinement halo model. We will therefore devote the rest of 
our discussion to an analysis of the diffusion halo models using the 
Y-ray data. These models appear not to be definitely ruled out by 
either the radio data or the age data (Bulanov, Syrovatskii and Dogiel 
» 1976, Ginzburg and Ptuskin 1976). 



IV. The Diffusion Calculation 

The model we shall employ to investigate the effect of a diffusion 
halo on the cosmic ray distribiition is often referred to as the 
flat diffusion model (Ptuskin 1974). In this model the cosmic-ray 
sources are assumed to be confined to a thin ('“sk of half thickness L^; 
the cosmic rays subsequently diffuse in a (possibly larger) disk of half I 
thickness L (L>L^). We assume cylindrical symmetry about the z axis 
which is perpendicular to the disk surface and the diffusing volume is 
taken to be of infinite extent in the radial direction r. The cosmic 
ray particles are free to diffuse in z and in r until they reach the 
surfaces z = + L at which point they escape freely; the appropriate 
boundary conditions for the particle density are therefore n(z=+L, r) = 0. 

In previous applications of this model (Ginzburg and Syrovatskii, 1971; 
Pacheco 1971, Ptuskin 1974, Bulanov et al. 1972, Guet and Pacheco 1973, 
Ginzburg and Ptuskin 1976 and refs.) the source is taken to be uniform 
"Lq j z _5 Lq and for all r. We, on the other hand are specifically 
interested in radial variatiotiS of the cosmic ray density since it is 
the radial dependence of the cosmic ray flux which, when folded with 
the observed distribution of interstellar gas, gives the observed 
galactic longitude variation of the cosmic y->^ay flux. 

In the following we shaM assume a constant scalar diffusion 
coefficient for the sake of simplicity. Different values of the 
diffusion coefficient in the z and r direction could be simply 
accommodated by using a different scale for the z and r distances but 
since the scale factors would be proportional to the square root of 
the diffusion coefficient values our results will not be too sensitive 
to this assumption. 
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The <1iffiision equation may be written in cylindrical coordinates 


as 

^ * -raT “ <') 

.where S(r,z) is the source strength density, D is the diffusion coefficient 
land we have made use of the assumed cylindrical symmetry. Since we will 
ultimately normalize the solution of (1) to the observed density of 
cosmic ray particles at the position of the solar system (r=19 kpc) the 
magnitudes of S and D will not enter and we may set 

S(r,z)/D = q(r) for jz 1 <_ Lq 
= 0 otherwise 

where various forms of q(r) will be tried and discussed later. 

In the appendix we show that the Green's function appropriate to 
(1 ) is given by 

G(z,z ; r,r') = - I cos(r^z) cos(k^z' ) I^(k^rJ K^(k^r^) (2) 

n=0 

where k^ = r^and r^ represent the smaller and the larger 

of the two variables r and r'. 

Since we are interested only in source functions that are uniform 

jin z between + we may integrate out the z' variable immediately 

to obtain for the cosmic ray density 
00 

n(z,r) = I cos(k z) f (r) (3) 

n=0 ^ 

where 

fn(r) = (Zn+I *) *^ 0 

= VV') ,,, 

+ Io(kn'') KQ(knr')]. 

r 
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The radial source functions, q{r), will be discussed in detail 
in the next section. They are mainly of the form 

q(r) = r^ exp (-Br) (5) 

where A and B have been determined by fitting to experimental data. 

Once a particular form for q(r) has been chosen. the integrals in 
(4) are performeo numerically to obtain the functions and 

the series in (3) is summed until a satisfactory convergence is 
obtained. 

As was previously stated the cosmic ray density function 
ncr(z»>^) obtained in this manner is normalized to the proper 
solar system value at z = 0, r = 10 kpc. This solar system value 
is obtained from the cosmic ray flux observed at earth by 
demodulation. 

V. The Cosmic Ray Source Distribution 

We consider various possible source distributions for cosmic 
rays ♦'O be connected with the formation of population I stars and 
their evolution to pulsars and supernovue. In this case, the 
distributions to be considered will be functions of galactocentric 
distance r (Burton 1976, Stecker 1976b). This distinguishes our 
diffusion halo models from those considered by Ptuskin (1974) and 
other workers at +fte Lebedev Institute (Ginzburg and Ptuskin 1976) who 
assume a uniform cosmic ray distribution throughout the galactic disk. 

We choose for examination the following radial distributions obtained 
from observational studies and analyses: 
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A) Supernova remnant distribution of Ilovaisky and Lequeux (1972). 

B) Supernova remnant distribution of Kodaira (1974). 

C) Pulsar distribution of Seiradakis (1976) (weighted and unweighted). 

D) Distribution of molecular clouds as given by Burton et al (1975). 

E) Distribution of ionized gas as given by Lockman (1976). 

In all cases, we consider the sources to be uniformly distributed 
in the z direction (perpendicular to galactic plane) over a disk of 
thickness 100 pc. Our results are insensitive to source width since, 
for any reasonable population I source distribution, the thickness of 
the source disk is much smaller than the halo thickness. 

The various experimentally determined radial distributions were fit 
with curves of the forn 

q(r) = r^ exp(-Br). 

The fit was a, standard least squares fit to the data points, and in all 
cases but one the points were given equal weighting. In the case of 
the pulsar distribution two fits were obtained: one with equal weighting 

and one where the points were weighted inversely as the square of the 
reported experimental error. As can be seen from the pulsar distribution 
and curves shown in figure 1 and from Table 1, including the experimental 
uncertainty in the fitting procedure can make a significant difference 
in the result. It should be noted that including the errors in 
fitting the pulsar data of Seiradakis (1976) yields a distribution 
that is almost identical to that obtained for the supernova remnants 
by Kodiara (1974). 
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is low in this central region and the path lengths are relatively 
small, the calculation of the Y-ray longitude distribution, as discussed 
in this section is not very sensitive to the exact value of q in the 
innermost regions of the Galaxy. 


VI. The Y-Ray Line Flux Longitude Distribution 

The cosmic ray specific intensity in the galactic gas disk is 
given by the solution of equation (1) as 

i| 

4 7T 


I6) 


It is convenient to express our length scales so that the distance 
between the sun and the galactic center (R^ ^ 10 kpc) is taken equal to 
1. The y-ray flux from Tf°-decay can then be expressed by 

I,(t) = i/ds q^(s) ( 7 ) 

where s is the distance along the line of sight and 

4^ nj(s)/dE 
j , k 

(see, e.g., Stecker 1971) where j and k represent the type of target 
and cosmic ray nuclei. 

For a cosmic mixture of H and He interstellar gas as target 
nuclei, and a value = I^^/I^.^ g , I^^ g being the demodulated flux 
of cosmic rays in the solar neighborhood (I « = I (r=l, z=0)) as given 
by Gloeckler and Jokipii (1967). Stecker (1970) gives a value for the 
Y-ray production rate in the solar neighborhood for energy greater than E. 

0 ,= ( 9 ) 

where qgj = (1.3 + 0.2) X 10"^^ 

(Stecker 1973, using the data of Comstock et al. (1972) gives an upper 
limit of qg^® ^ 51 ^ ig-25 variables defined in Figure 2, 
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and the assumed distribution c(r,e,z) = c(r,z=0) = i{r)t 


n,^(r,0,z) = n^(r,z=0) = n^(r) 
equation (7) becomes 


I U) = 

^ 4ti 


/ ds n^(r) C(r) 


( 10 ) 


where r = (s^ - 2s cos-^ + 1)® 

In this particular case, we wish to calculate the rray flux above 
100 MeV as would be seen by the SAS-2 detector integrated over 
galactic latitudes - 10° b ^ 10°. We consider the hydrogen gas 
to be made up of two disks, an atoi.’ic hydrogen disk with a half- 


width h^ =1.1 X 10"'^ Rq (110 pc) and a molecular hydrogen disk with 

h2 = 5 X 10"^ Rg (see, e.g. Burton 1976). In the outer galaxy 

(1 <_ r ^ 1.5), this disk width is larger by a correction factor 

w(r) = 1.8r - 0.8 (Baker and Burton 1975) which we take account of by 

using w(r) as a weighting factor in the integral (with w(r) = 1 for 0jr<^l. 

With a mean atomic hydrogen density n^j(r) as given by Burton et al . (1975) and 

mean molecular density n^ (r) as given by two determinations: (a) 

H2 

from Burton et al. (1976) and(b) from Scoville and Solomon (1975) 
normalized by Stecker et al. (1975). We then obtain the formula 
for the line intensity 


>100 

-r 

I U) = 


>100 




10° 


hjcot b 

i2COt b 

db 

0° 

-I 

ds n^j(r)c(r)w(r)+?’ 
. 

) 

ds n^^(r)f;(r)w(r) 

0 


( 11 ) 


where qg j is the total y-ray production rate per !< =itom above 100 MeV, 


>100 >100 


’9,T 


= q 


>100 


>100 


0TT 


^0,Bremstrahlung ^0, Compton 


( 1 ?) 


== 1.33 X 10"^^ 5“1 (Stecker 1977) 
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(Strictly speaking, the Compton y-r&y production rate does not scale 
with but this contribution is only about 5% of the total production 
rate so that only a small error is introduced.) 

The theoretical fluxes calculated using equation (11) for 

varius diffusion models and source functions were compared with the most 

recent SAS-2 data (Kniffen, Fichtel and Thompson 1977) with point 

2 

sources subtracted out. In each case, a x test for goodness of fit 
was made and the statistical probability of obtaining the SAS-2 longitude 
distribution was obtained, assuming the model represented the real y-ray 
longitude distribution. The tests were made over the longitude range 
10° L 90°, since detailed CO surveys of the molecular cloud 
distribution are only available in the range 0° < £ < 90° and the range 
0° <_ 10° is expected to contain a significant contribution from 

Compton produced y-rays which is not included here (see Stecker 1977 
and Section VIII). indeed a"*! smooth models of the large scale galactic 
distribution of y-rays do not reproduce features in the observed longitude 
distributions for I > 120° which must be produced by nearby fluctuations. 
This is supported in the SAS-2 contour maps for the anticenter region 
(Hartman et al. 1976). 
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VII. Goodness of Fit of the Theoretical Longitude Distributions - x 


After the longitude distributions of ■rr'ay flux were calculated 

as described in the previous section they were compared with the 

fluxes observed by SAS-2 (Kniffen et al . 1977) to see which, if any, 

gave reasonable fits to the data. We determined the goodness of 

2 

fit by use of the x statistic of K. Pearson (Cramer, 1946). The 
advantage of the x statistic over a simple sum of squares of the 
residuals is that it enables an absolute measure of goodness-of- 
fit to be given for a particular theoretical distribution. 

If the data are broken up into R bins of width 2.5 degrees 
of longitude, theory gives the expected number of photons observed 
in that bin to be = a^F^. where F^. is the predicted flux in that 
bin and a^. is an instrumental factor made up of such things as 
detector sensitivity, look time, etc. The SAS-2 data, on the 
other hand, are reported in terms of a flux f^ so we may assume that 
the observed number of photons was = a^f^ with the same set of 
The reported errors are primarily counting errors so that 

We may eliminate the to obtain 


^ 1 

6N 

z 

' “i 


oi 

N . 


ffx.) 

01 




^ni = 

“i " 

f 

1 ( 

We 

may 

now I 

s' 

R 

= Z 



- / N 


01 


R 


c=l 


’ei 


.i. fh (^^1 


c=l cTfFf 
2 


' fi 


Test 


■s. 


(13) 


(14) 


Where the final form is written soley in terms of the theoretical 
fluxes, the observed fluxes, and the experimental errors for the 
fluxes. 

If the observed counts are normally distributed about the 

‘ 2 

"true" expected value then it can be shown that the quantity s 

of (3) is distributed with a x -distribution with R degrees of 

freedom (Cramer, 1946). In fact the have a Poisson distribution 

but if the mean value is large enough the Poisson distribution 

is well approximated by a normal distribution so one would 
2 2 

expect the sum s to have approximately the x -distribution. From 
2 

tables of the x -distribution one can read the probability that the sum 
given in (3) would have a value greater than or equal to any given value. One 
can test the hypothesis that a given theoretical model represents the true 

expected flux values by noting the orobability that one would obtain a value 
2 

of X as large as the one actually obtained upon comparison of the model with 
the SAS-2 data. If this probability is too small e.g. <5% then 
it is usual to assume that the hypothesis is unlikely to be correct 
and the theoretical curve is rejected. 

» Whereas a theoretical model that yields a higher probability than 

an alternative model is usually to be preferred, one must apply this rule 

with caution. If one finds for example that there is a 99% chance 
2 

of having a x as big or bigger than the one obtained, this 

indicates a very close fit of the data to the theoretical model; yet 

it also indicates that there was only a 1% chance of obtaining this 

close a fit to the "correct" model. This is in itself quite suspect. 

The general rule to be applied in interpreting these probabilities 

il that any values near 50% should be considered a good fit, for this 

is what one would expect most of the time with the "true" model. 
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Probabilities sensibly smaller than this (i.e. < 10%) should be 
considered as poorer fits and probabilities less than 1% are quite 
unlikely and indicate that the model being tested is probably 
not correct. 

In the next section we shall present the probability values obtained 
for the various models we h^e considered and discuss what we believe 
are the implications as to the existence and extent of a significant 
cosmic-ray halo. 

VIII Results and Discussion 

Tables 2 and' 3 show the statistical probability that the y-ray 

data will be produced with the observed intensities and errors in 

the SAS-2 data in the 10° <t <90° range assuming that each 

theoretical model gives the "true" longitude distribution. These 

2 

probabilities are obtained from the x test. The models are specified 

by the source functions described in Section V and the L values 

given in the tables. In the one dimensional approximation (diffusion 

^ 2 2 

perpendicular to an infinite plane (d n/dz =0, n^^(l) = 0), 

Icr(z) = (z<L) (15) 

Thus, the "size" of each halo model is smaller than L if one speaks 
in terms ofal/ehalf thickness (0.63L) or a halfjwidth half maximum 
(HWHM 0.5L). 
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TABLE 2. SAS-2 DATA PROBABILITY OF FIT {%) TO DIFFUSION 

HALO MODELS* 


Source 


Model 

No Halo (L=0) L=5kpc 

L=10 kpc 

L=20 kpc 

SN (Ilovaisky and 
Lequeux 1972) 

32;16 

3;0.06 

« 0.1 

<<0.1 

SN (Kodaira 1974, 
Broken Linear 
Fit) 

45;67 

16;3 

« 0.1 

« 0.1 

SN (Kodaira 1974) 

42; 39 

6;0.25 

<< 0.1 

<< 0.1 

Pulsar (Weighted, 
Seiradakis 1976) 

44; 47 

7.3;4 

• 

A 

A 

O 

« 0.1 

Pulsar (Unweighted) 

« 0.1 

'v.0.1;22 

45; 69 

30;12 

CO Cloud Distribution 
(Burton et al . 1976) 

<< 0.1 

1;38 

48;62 

25;8 

Ionized Gas Distribution 
(Lockman 1976) 

«0.1 

« 0.1 

3; 50 

51;73 

*First number is obtained 

using N^^ 

from Burton et al . 

(1976). Second 

number 


is obtained using from Scoville anj Solomon (1975) as normalized by Stecker 
et al. (1975), E.g. B;S. 

TABLE 3 

PROBABILITIES OF FIT (%) USING THE KODAIRA SOURCE FUNCTION 


L 

0 

1 kpc 

3 kpc 

5 kpc 

10 kpc 

with Burton et al. n^ 

”2 

42 

40 

19 

6 

<< 0.1 

with Scoville and Solomon 

39 

31 

4 

0.25 

<< 0. 1 


"H2 
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It can be seen from the tables that if one assumes a source 
distribution as given by the supernova remnant or pulsar (weighted) 
distribution, fits with probabilities greater than ^S% are found 
only for models with small halos (L< 5 kpc or, from (15), a "scale height" 

< 3 kpc) with the most probable models being the ones with very thin halos 
or no halos (L=0). This is also consistent with the radio data (Baldwin 
1976) and some isotope measurements (Hagen et al. 1976), although the radio 
and isotope measurements do not in themselves rule out extensive diffusion 

. halos (Ginzburg and Ptuskin 1976; see also section II). 

If one wishes to save the extensive halo hypothesis, it appears to 
be necessary to relate the source of cosmic rays to a sharper distri- 
bution such as that of molecular clouds or ionized gas rather than the 
presently favored supernova-pulsar distribution. One such hypothesis 
which has been discussed is that of large scale first order Fermi 
acceleration associated with the density wave theory of star formation 
(Stecker, et al 1974, Puget and Stecker 1974). 

It can also be seen from the tables that it makes no qualitative 
difference to the results whether one adopts the molecular cloud distri- 
bution given by Burton et al (1976) or that given by Scoville and 
Solomon (1975) in computing the y-ray fluxes. 

Finally we have also examined with the x test, some galactic 
Y-ray models discussed in the literature. The model considered by 
Stecker (1975) is found to have a 39% probability in the longitude range 
10°<^£ ^90°. The more elaborate model of Kniffen, Fichtel and Thompson (1977) 
based on spiral arm enhancements, leads to only a 5% probability, 
reflecting a lack of detailed correlation between the arm features 
of the model and the peak in the observed data in the 10° <_£ <_90° range. 
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Consideration of other longitude ranges may show a better 
correlation, however. Extrapolating from models which we have investigated 
with large L, the constant cosmic ray model of Fuchs, Schlickeiser and 
Thielheim (1976) would provide a poor fit (probability <<1%) to the data. 
This is expected as it would correspond to the case L-*^. 

Figure 3 shows the cosmic ray distribution in the plane of galaxy 
obtained using the broken linear fit to the Kodaira SN distribution 
for diffusion halos with various values of L (all normalized to 1 at 
10 kpc (r=l)). Figure 4 shows contours of constant cosmic ray 
intensity in r and z obtained from this model for L=l,3 and L=10 kpc. 

The L=1 and 3 kpc cases can fit the v-ray data, but the L=10 koc case is 
ruled out by the y-ray data. 

Figures 5 and 6 show the y~ray longitude distributions calcu- 
lated using the weighted pulsar distribution for a source function 
with no halo (L=0) and a 10 kpc halo. The Scoville and Solomon 
n^ was used in these models. They are compared with the latest 
SAS-2 longitude distribution with point sources removed (Kniffen, Fichtel 
and Thompson 1977). The L=0 case has a probability of 47%. The L=10 Kpc 
case has a probability of fit of 'v- 10"®. 

Figure 7 shows the model of Figure 5 with an additional contribution 
from Compton interactions as calculated from the model of Stecker (1977) 
with he = 500 pc. The value of he used was obtained by a minimum 
fit, normalizing the theoretical Compton longitude profile. In this no- 
tation, h would correspond to the half-thickness of the electron disk 

(L ) and £ = I /I« is the normalized cosmic ray electron flux in the 
e e 0 

galactic center region, i.e., if e = 1* = 0.5 kpc at the galactic 

center. 
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One may attempt to determine the thickness of the cosmic ray electron 
halo directly from the latitude observations of galactic ^-rays, as has 
been attempted by Schlickeiser and Thielheim (1977). These ithors have 
estimated a value for of 3.2 kpc. They obtained this value by taking 
the SAS-2 data, subtracting out an estimated extragalactic component and 
an estimated galactic component and then assuming that the remainder is due 
to Compton interactions of cosmic-ray halo radiation with the blackbody 
and starlight radiation fields. However, it is not clear how these authors 
obtained a value for the extragalactic component and separated out a halo 
component. Unfortunately there are also large uncertainties in all of the 
estimates involved using this method. Stecker (1977) has pointed out that 
the SAS-2 data at high galactic latitudes is consistent with the sum of 
galactic and extragalactic contributions with no significant halo contribu- 
tion. However, some halo contribution is certainly also consistent with 
the SAS-2 data, so that we examine ^his question in more detail here. 

We assume the SAS-2 high latitude flux to be made up of three com- 
ponents, a flux from the galactic disk, a Compton flux from a thick disk 
or halo of half-thickness L^, and an extragalactic component: 

^SAS " ^gal ^ ^C ^ex (16) 

In our case, we take the extragalactic flux to be determined by an 

_3 

E power-law fit in accordance with cosmological pion-decay theories of 
the diffuse background (Stecker 1975b). We then obtain, for y-ray energies 
above 100 MeV. 

Isas = ^.9 + 0 . 4 ) x 10 "^ (Fichtel et al. 1975) 

Igai 0-6 X 10 ^ cm'^s'^sr’^ (Stecker 1977) 

Igj^ = 0.8 X 10"® (Stecker 1977) 
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Thus, by combining (16) and (17), one obtai :s a value for 1^, of 

fi -? 1 -1 

(5+4) X 10 cm s sr . The uncertainty given here is a lower limit, 
not reflecting the uncertainty in the theoretical estimates but only the 
uncertainty in the data. The value for is then given by 

“■e ' ^00 'c ^ 

‘'e.c 


( 18 ) 


where b is the mean latitude at which the SAS-2 flux was observed. 

The Compton production rate for the solar vicinity for energies 
greater than 100 MeV, » is given from recent estimates as 

6 X 10'^^ (Stecker 1976a) 

■ 5.5 X lO’^^cm^s'^ (Piccinotti and Bignami 1976) (19) 

‘ 8.5 X 10"^^ (Stecker 1977) 

The uncertainties depend on various estimates of the optical and infra- 
red photon fields and on the exact sh?pe of the electron spectrum. We 
adopt here a value of (7 1.5) x 10'^^ cm'^s*'\ again probably under- 

estimating the uncartainty. We also take for a mean latitude, b = 40°. 

We then obtain from equation (18) an estimate for the mean half-thickness of 
the electron halo of 


Lg = 2 + 2 kpc. 

The radio measurements and the cosmic ray age measurements also 
tend to favor, or at least be consistent with a thick disk or thin halo 
with L = (1-3) kpc (Ilovaisky and Lequeux 1972, Baldwin 1976, Stecker 1977, 
see also section I). Our present results on the cosmic ray nucleon halo, 
using a supernova-pulsar type source distribution are also pe*'fectly con- 
sistent with this thick disk model, as can be seen from Table 3. 
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One can use the values obtained for L from the diffusion model with 
a supernova-pulsar ^ype source distribution, i.e. L<5 kpc in order to 
obtain values for the diffusion coefficient D and the mean-free path for 
diffusion x. These quantities can be evaluated from the relations 


D = L^/2t 
X = 3 D/c 

where x is the cosmic ray age (Ginzburg and Syrovatskii 1964). 

Estimates of D and x obtained from equation (20) are shown 
4 and 5. 


( 20 ) 


in Tables 


Table 4. Values for D (cm^s"^ ) 
x(10®yr) 


L(kpc) 

3 

10 

20 

1 

4.5 X 10^® 

1.5 X 10^® 

7.5 X 10^^ 

2 

1.8 X 10^^ 

6.0 X 10^® 

3.0 X 10^® 

4 

4.1 X 10^^ 

1.4 X 10^^ 

6.8 X 10^® 

5 

1.1 X 10^° 

3.8 X 10^^ 

1.9 X 10^^ 


Table 

5. Values for X (pc) 




t( 10^ yr) 


L(kpc) 

3 

10 

20 

1 

1.5 

0.5 

0.25 

2 

6.0 

2.0 

1.0 

3 

13.5 

4.5 

2.25 

5 

37.5 

12.5 

6.25 

Val ues 

of L greater than 

5 kpc are considered 

to be highly 


unlikely (see Table 3). 
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APPENDIX 


Consider the eigenvalue equation 


iV . 1 

3 r 


3 n 
3 r 


- An = 0 


(Al) 


The above equation is separable in z and r and the solutions 
that fit the appropriate boundary condictions (n=0 at z = + L, 
integrable in r) are 

n^(z,r) = cos(k^z) JQ(k.r) 
where k^ = 
and A=-(k^"+k ) 

Since the diffusion operator is self adjoint, the eigensolutions form 
a complete set and it is straightforward to verify that 


^ cos(k^z) cos(k^z' ) = ^ (6(z-z ) + 6(z+z' )) (A2) 

(Since we are interested only in models that are symmetric 
in z, we have employed only the symmetric eigenfunctions resulting in 
the symmecric delta function.) 

Ob 

dk k J (kr) J (kr‘ ) = ^ 6(r-r' ) (A3) 

^ 0 0 r 

If one has a complete set of orthonormal eigensolutions of the 
) equation 

L Uj • XUj (A4) 

Such that 
» 

Ua(x) Ua(X') = 6(x-x') (A5) 

4 

A 
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1 

i 

I 

Then the Green's function defined by 

L G(x-x' ) = 6{x-x' ) , 


is constructed from the eigensolutions. 



< (Cushing 1975, Jones 1970). 


For our problem we thus have 

00 

G(z,2',r,r') = - r I cos(k.z) cos(k_z')f dk k Jn(kr)Jo(kr' ) 

'• n=0 " "I k?+k„2 (A7) 

The integral over k may be performed to yield 

G(z,z', r,r') = - cos(k^z) cos(k^z' )lQ(k^r ) K^(k^r^) (A8) 

where r< (r>) is the lessor (greater) of r and r' and Iq and are 
the modified Pessel functions of the first and second kinds 
respectively. 
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Figure Captions 


Figure 1 



Figure 2 

» 

Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


The pulsar distribution in the Galaxy obtained by Seiradakis 
(1976) together with best fits to a distribution function of the 
form r^exp(-Br) both unweighted (curve A) and weighted according 
to the error bars as shown (curve B). 

Definition of distance variables used in the y-ray flux 
equations. (G.C. = galactic center) 

Cosmic ray distribution in the galactic plane (z=0) using a 
broken linear fit to the supernova remnant distribution of 
Kodaira (1974) and various diffusion halo sizes L (in units of 
Rq = 10 kpc). 

Contours of constant cosmic ray intensity in the r-z plane for 
the source model of Figure 3 with L=l,3 and 10 kpc. 

The Y-ray longitude distribution obtained using the weighted 
pulsar source model with L=0 (no halo source dominated model) 
compared with the SAS-2 data. (Probability of fit equals 47%). 

The distribution of '^^s obtained from Scoville and Solomon 
(1975) and Stecker et al. (1975). Using the Burton et al. (1975) 
distribution for n^^ or the supernova remnant source distribution 
of Kodaira (1974) also gives a good fit (see Table 2). 

The Y->"ay longitude distribution obtained using the weighted 
pulsar source model with an L=10 kpc diffusion halo compared 
with the SAS-2 data. (Probability of fit 10"^). 

The Y-ray longitude distribution obtained using the same model as in 
Figure 5, but with an additional Compton Y-ray flux as discussed 
in the text. The SAS-2 data are again shown for comparison. 
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